Interaction-induced adiabatic non-linear transport 
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We calculate the time-dependent non-linear transport current through an interacting quantum dot 
in the single-electron tunneling regime (SET) . We show that an additional dc current is generated by 
the electron-electron interaction by adiabatic out-of-phase modulation of the gate and bias voltage. 
This current can arise only when two SET resonance conditions are simultaneously satisfied. We 
propose an adiabatic transport spectroscopy where lock-in measurement of a "time-averaged stability 
diagram" probes interactions, tunnel asymmetries and changes in the ground state spin-degeneracy. 
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Introduction. Transport through nano-scale devices 
modulated by time-dependent externally applied electric 
fields is an active field of research important for trans- 
port spectroscopy and manipulation of the charge and 
spin degrees of freedom in nano-structures, see, e.g., 
A particularly gentle way of time-dependently probing 
a system is through "adiabatic pumping" [2Hj]. Here a 
finite dc current is generated in the absence of an ap- 
plied bias by a weak, low frequency periodic modulation 
of system parameters. Adiabaticity in a transport situ- 
ation means that many electrons visit the system dur- 
ing one cycle of the driving with frequency f2 and that 
the modulation is too weak and too slow to excite the 
device by direct state-to-state transitions. The current 
generated this way crucially depends on which subset of 
parameters is modulated, on the working point about 
which the modulation takes place and on interactions, 
which are of special importance in nano-scale devices. 
For strictly adiabatic pumping one needs to vary at least 
two parameters, single parameter pumping requiring a 
higher frequency [f|. Among the various combinations 
of parameters studied so far, the modulation of the ap- 
plied bias has received little attention Most works 
have considered small deviations around an equilibrium 
working point where no steady state current is flowing. 
Adiabatic modulation around a non-equilibrium trans- 
port state induced by a static non-linear bias voltage has 
been explored only for systems with negligible Coulomb 
interaction [4J, motivated by experiments with surface 
acoustic waves 0j- Non- linear bias voltage and Coulomb 
interaction have received little theoretical attention in 
the adiabatic regime. Limited to an equilibrium work- 
ing point, some works have studied interacting quantum 
dots [1, Q and wires . Including the effect of strong 
interactions beyond the mean-field picture is a challenge, 
since the powerful scattering matrix approach Q breaks 
down here. Generally, one expects the additional non- 
equilibrium introduced by a static dc-bias voltage, in 



combination with strong electron-electron interactions, 
to strongly modify the pumping, providing novel oppor- 
tunities to investigate and control transport properties of 
nano-scale devices. 

In this letter we propose a new scheme for transport 
spectroscopy of interacting systems using adiabatically 
time-dependent electric fields. We analyze an interact- 
ing quantum dot in the SET regime, adiabatically driven 
by out-of-phase gate and bias potentials. In contrast 
to previous works, the applied bias can be arbitrary, 
i.e., we modulate the parameters around a steady non- 
equilibrium state supporting a finite dc current. We 
show that the strong local interaction generates an ad- 
ditional adiabatic dc current, which is identically zero 
without interaction for any value of the applied volt- 
ages and magnetic field. We propose to use this effect 
as a tool for non-linear transport spectroscopy which can 
be measured using lock-in techniques. The adiabatic dc 
current is non-zero only when two conditions for single- 
electron tunneling are simultaneously satisfied. Plotted 
as function of the time-averaged gate and bias voltage, it 
gives rise to a new type of "stability diagram" . Further- 
more, we show that in an external magnetic field lifting 
the spin-degeneracy, the adiabatic modulation only gives 
rise to transport effects in the regime of non-linear bias, 
which qualitatively distinguish between different junction 
asymmetries. 

Model. We consider a quantum dot weakly coupled to 
two electrodes as sketched in Fig. [TJa). The gate and 
bias voltage are modulated with frequency Q around the 
working point specified by the voltages V g and Vt>: 

V x (t) = V x +SV x sm(nt + ip x ), x = b,g. (1) 

We consider the important case where a single orbital 
level with strong Coulomb interaction U is relevant for 
transport. We denote the spin-resolved dot number op- 
erator by n a = d' a d a , where the spin a =t, 4- is quan- 
tized along the external magnetic field (if present). The 



2 



Hamiltonian reads H^>(t) = e a {t)n a + Un^n^. The 
energy of an electron created by d\ equals e a (t) — 
— aV g (t) + aB/2 using the shorthand a = ±1 for spin 
t,4- Importantly, the time-dependent gate voltage V s (t) 
capacitively modulates this energy with lever arm a < 1. 
Furthermore, B is the Zeeman energy in units e = H = 
fcs = 1. The many-body eigenstates of Ho(t) are |0), 



is the electron temperature, the result reduces to Fermi's 
Golden Rule. In Eq. (j3]) the voltages arc replaced by their 
instantaneous values at measuring time, V x — > V x (t), 
x = b, g, resulting in a parametric time-dependence indi- 
cated by the subscript t. Hence, would be the time- 
dependent steady state of the dot if the system was able 
to follow the parameter modulation instantaneously. By 



\a) = d+|0) with a =|,| and |2) = d\d\\$) with energies inserting p(t) = pj 1 ' +p| a) in Eq. ©, using Eq. ©, and 
0, e a (t), J2a £ cr(^) + U, respectively. The time-dependent 
bias Vb{t) enters through the electro-chemical potentials 
fi r (t) = ±Vb(t)/2 of electrodes r = L,R, which are de- 



doing a systematic lowest order expansion in Q/T <C 1, 

fa) 

we hnd for the hrst correction p). to the instantaneous 
(0. 



scribed by H E (t) = J2r,k,a ( e k + A*r(*)) c\ ka c rka . Finally, 
Ht = J2 r k a trCrkad\ + H-C- describes the tunneling be- 
tween the dot and the electrodes, with tunnel coupling 
strength T r = 2irp r \t r \ 2 , where t r is the amplitude and 
p r the density of states of lead r = L, R. We define 
r = ^ r T r . We consider here the important case where 
the transport is affected dominantly through the modula- 
tion of the energy level positions s & (t) and the bias energy 
window fii,(t) — /iR,(t). The change in the tunnel coupling 
is of negligible importance. This is the typical situation 
in non-linear transport spectroscopy of quantum dots. It 
holds in particular for small amplitude modulation of the 
voltages considered here. 

Retarded occupations and transport current. The to- 
tal Hamiltonian H(t) = H E (t) + H B (t) + H T thus 
contains strong interaction on the dot, an adiabatic 
time-dependence and non-equilibrium introduced by the 
non-linear bias voltage. Within the framework of the 
real-time transport theory the time-dependent occupa- 
tion probabilities of the many-body dot states, p(t) = 
(po(t),pi(t),pf(t),p2(t)) can be shown to obey the ki- 
netic equation [9( 



reference solution p t : 
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dt'W(t,t')p(t') 



(2) 



The kernel, W(t, i'), accounts for changes of the dot oc- 
cupations due to electron tunnel processes to/from the 
electrodes. Although, it explicitly depends on both time 
arguments t and t' (in contrast to the time-independent 
case) it can be calculated perturbatively for slowly vary- 
ing fields 0] . Here we restrict ourselves to the lowest or- 
der contributions in both the tunneling coupling (single- 
electron tunneling (SET)) and in the time-dependent 
perturbation of external system parameters (adiabatic 
driving). We consistently solve the kinetic equation by 
expanding it around the instantaneous reference solution, 
Pt , defined by 



= W t (i) p«. 



(3) 



The instantaneous kernel and its zero-frequency Laplace 
transform W 4 (i) = lim^o dt'W^(t - i')e*rt*-0 are 
evaluated using stationary transport theory in lowest or- 
der in r. In this limit of weak coupling, r -C T, where T 
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This correction accounts for the actual delay suffered by 
the system due to the finite rate of sweeping the volt- 
ages. Further corrections to this adiabatic approxima- 
tion can be neglected if in addition a5V g ,5Vb <C T 2 /il. 
The time-dependent steady state including the retarda- 
tion is uniquely determined by equations ([3]) and (j4|. 
together with the normalization conditions e T pj I ' ) = I 
and e T p< a) = with e T = (1, 1, 1, 1). 

The time-dependent current flowing from lead r = L, R 
into the dot is found in a similar way Q and can 
be decomposed into two corresponding parts, I® and 

l[ a } . Here 1$ is the adiabatic correction to the cur- 
rent due to the retardation of the system, i.e., it van- 
ishes in the limit 5V£l — » 0. The central quantities dis- 
cussed in this paper are obtained when averaging the two 
current contributions over one entire modulation cycle 
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J^r J Q n dt I t fj. . Here, If' equals the dc cur- 
rent one would measure for time-independent voltages 
equal to V s and Vb- Plotting dl^/dVb as function of 
these voltages, one obtains the standard Coulomb block- 
ade stability diagram [l[ . The quantity of central interest 

-fa) 

here, J-r , is the additional dc current component due to 
the retardation of the quantum dot state. This quan- 
tity can be obtained, e.g., by subtracting from the total 
measured time-averaged current its zero frequency limit. 

For the time-dependent adiabatic current, we obtain a 
central analytic result 



r (a) 



(r r + 7r )(r-7) + /3 r ff d_ 

T 2 - 7 2 + /3 2 dt 



n) 



o (r r + 7r )/?-&.(r + 7 ) d (i) 



r 2 - 7 2 + p 2 dt 

The average instantaneous charge, (n) 



(5) 



(i) 



~^2t ai a l^)Pt \ are found from 



m% and spin, (S z )f> 
Eq. ([3]). Although the time-dependent adiabatic cur- 
rents depend on the junction r where they are evaluated, 
their time averages are related by charge conservation, 
J2 r %r = 3?( n )* l) ' gi vm g J2 r ^ = 0. The prefactors in 



3 



Eq. ((5j) contain 

lr (t) = ir r J2 [/ (*))-/ (Car (* ) + U)] (6) 

(T 

/3 r (t) = ir r £ a [/ (e CTr (i)) - / (e ar (t) + U)] (7) 

a 

and their sums by 7 = $^ r 7 r , an< ^ P = Er^n where 
(T = ± (corresponding to t, 4-) and r = L,R. All these 
quantities depend on time through the distance to res- 
onance e ar (t) = e a (t) — /j. r (t) in the arguments of the 
Fermi-function f(tu) = (exp(cj/T) + l)" 1 . From Eq. © 
we infer a necessary condition for a non-vanishing time- 
averaged adiabatic current which also holds for more 
complex systems: 7^ can only be non-zero if two SET 
resonance conditions are satisfied simultaneously. If only 
a single resonance condition is satisfied (effectively this is 
single-parameter pumping), (|5|) is a total time-derivative 
of a periodic function, resulting in a zero time-average. 
The resonances in i£ are thus located at resonance line 

crossings of the standard dl^'/dVb map. 

Interaction-induced dc current. The central result 
of the paper relates to the prefactors in Eq. 
Since the tunnel rates r r and T = ^ r r r are time- 
independent it is clear that the adiabatic dc current 
is generated by the Coulomb interaction U. Indeed, 
since 7r = fi r = for U = the adiabatic cur- 
rent 1$ = (r r /r) (d(n)l i) /dt + 2d(S z }£ ) /dt) is a total 
time-derivative, which, integrated over a period, yields 
Ir^ = 0. We emphasize that in this case the current Ir^ 
vanishes identically for any value of the time-averaged ex- 
ternal voltages and of the time-independent tunnel cou- 
plings and external magnetic field. We now discuss the 
voltage dependence of the adiabatic dc current in the 
experimentally important regime of strong local interac- 
tion U >T>T. We fix the direction of the modulation 
cycle by taking tp\> = and ip g = —n/2, for which the 
adiabatic dc current is maximal, and time-average the 
current numerically. We first focus on the case of zero 
magnetic field for which e a (t) = e(t) = —aV g (t) is inde- 
pendent of spin a. Therefore (3 r = and Eq. ([5]) simpli- 
fies to 7 t (a) = (r r + 7r )/(T + 7) d{n)f/dt. In Fig. O^b) 

we show a time-averaged stability diagram, i.e., i£ plot- 
ted as function of the time-averaged gate and bias volt- 
age. In contrast to the standard (dz£ /dVb) stability 
diagram in the inset of Fig. |Tfb) this map of pumped 
current indeed shows resonant enhancements only at dis- 
crete points of size oc T where two SET resonances 
meet. Most prominent are the two charge degeneracy 
points (aVg,Vb) ~ (0,0) and (U, 0) at which the adia- 
batic dc current has opposite sign and maximum am- 
plitude. We now explain the microscopic origin of the 
positive sign of the adiabatic dc current at the degener- 
acy point (aVg,Vb) ~ (0,0) for symmetric tunnel cou- 
pling Tl = Tr. For t <E (0,7r/f2) the adiabatic current 




FIG. 1. (Color online), (a) Sketch of the transport setup, (b) 
Adiabatic dc current through the left junction, 7^ , as func- 
tion of the time-averaged gate and bias voltage for Coulomb 
interaction U = 30T, zero magnetic field and F — 0.5T, A = 
0.25 and driving parameters Q = 0.1T, a8V g = 5V b = 0.5T. 
We plot 7^ scaled to the maximal absolute value at the de- 
generacy points, = ff ^^^^^r, consisting of the 
frequency, a coupling asymmetry factor, and the ratio of the 
voltage phase-space factors to the thermal energy window. 
Inset: differential conductance, dl^ /dVb, versus V s and Vb. 

(c) Voltage modulation cycle around the degeneracy point 
(aV g ,Vb) ~ (0,0) with dashed resonance lines e(t) = fi r (t). 

(d) Same as (b) for finite applied magnetic field, B — 10T. 

through the left junction is positive, 7 t ^ > 0, whereas 
in the second half of the cycle 7^ < 0. This is be- 
cause r + 7 « r(l + J2 r /( e (*) - M*))/2) and d(n)f ] /dt 
are symmetric and antisymmetric functions of the time 
t. The time-average 7^ is nevertheless non-zero due 
to the factor T L + 7L w T/2 + Tf(e(t) - yU L (i))/2 in 
the numerator. Clearly, since the first term is constant, 
the non-zero time-average comes from the contribution 
cx f{s(t) — fJ*L.{t)) which is non-zero for times for which 
e(t) < /J,h(t) (red part of cycle in Fig. [He)). One thus 
samples predominantly the loading parts of the cycle 
where d(n)[^/dt > (shaded in Fig. [He)), where an 
excess of electrons tunnels onto the dot through the left 
junction. Therefore the adiabatic dc current is positive. 
Similarly, one finds for the point (aV s , Vb) ~ (U,0) the 
opposite adiabatic dc current due to the negative sign 
of the second term in Eq. For asymmetric rates 

Tl 7^ Tr the time-dependence of T + 7 becomes impor- 
tant as well, but does not alter the sign of the adiabatic 
dc current. In a magnetic field B ^> T the adiabatic dc 
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FIG. 2. (Color online). The Iff along the line e = /^r 
(dashed line in Fig. Hid)), as function of the tunnel coupling 
asymmetry parameter A, all other parameters being the same 
as in Fig. O^d) and L is taken for A = (r L = Tr). 



current plotted in Fig. QJd) is completely suppressed in 
the linear response regime Vb <C T. Indeed, in this limit, 
j3 r and l[ a } has zero average (e.g. around aV g <C T 



it is I, 



(a) 



(r r /r) dp^ ydt). In general, the interaction 



breaks the symmetry of loading and unloading parts of 
the cycle. The magnetic field, however, restores this sym- 
metry in the linear response regime by lifting the spin de- 
generacy. Therefore, Ir is suppressed for Vb <C T even 
though U 7^ 0. This is to be contrasted with the stan- 
dard dI^/dV h map shown in the inset, where in linear 
response the conductance shows the Coulomb oscillation 
peaks. Only at a finite voltage Vb = B, where the spin- 
excited state becomes available, the loading-unloading 
symmetry is broken again and is restored. 

Coupling asymmetry. Asymmetric rates induce addi- 
tional features in Fig. (Ub) and Fig. Did) a * finite bias 
Vb| = U and U — B, respectively. In Fig. [2] we plot 
along the dashed line in Fig. HJa) as function of the 
coupling asymmetry, quantified by A = (IY — Tr) /P. 
Strikingly, the two high bias features marked (3) and (4) 
are qualitatively sensitive to the coupling asymmetry: if, 
e.g., resonance (3) is negative (positive), then Tl > Tr 
(Fl < Tr). Quantitatively, for A > the adiabatic dc 
current resonances marked (2) and (3) deviate from the 
"bare" resonance positions (A = 0) by a shift which de- 
pends linearly on the temperature T [111 . One can thus 
sensitively probe the coupling asymmetry. 

Adiabatic spectroscopy. Our results generalize to quan- 
tum dots with more complicated states and spectra: 
without interaction, the adiabatic dc current vanishes 
in leading order in T and fl. Therefore, measurement 



eral, the occurrence and sign of adiabatic dc current at 
a charge degeneracy point can be tied to the change in 
spin-degeneracy in the ground state: the sign is posi- 
tive (negative) if the ground state spin-degeneracy in- 
creases (decreases) with the quantum dot charge and it 
vanishes if there is no change. The time-averaged stabil- 
ity diagram thus directly reveals non-degenerate ground 

~( a) 

states if 7 L vanishes in the linear response regime. This 
may be interesting, e.g., for transport through magnetic 
molecules with high spin degeneracies or in carbon nan- 
otubes where both spin- and orbital-degeneracies play a 
role. Another important aspect of the proposed spec- 
troscopy is that the effects of "spurious" modulation of 
the barrier can be clearly identified experimentally. As 
shown in Ref. Q , a modulation of the gate voltage and of 
the barrier (instead of the bias voltage) results in an adia- 
batic dc current which is symmetric with respect to rever- 
sal of the time-averaged gate voltage aV g — > U — ceV g , in 
contrast to the antisymmetric shape found here. The pro- 
posed spectroscopy does furthermore not rely on quan- 
tum fluctuation effects and can therefore be observed 
readily in weakly coupled devices at moderate tempera- 
ture and low driving frequency. We have checked that 
the corrections from next-to-leading order tunnel pro- 
cesses (r 2 ) to the effects discussed here are quantitative 
and small, even for T ~ T. Importantly, even when in- 
cluding these corrections the adiabatic dc current still 
vanishes exactly for zero interaction. By measuring the 
proposed time- averaged stability diagram one thus gen- 
tly probes junction asymmetries and strong interaction 
effects. This may prove valuable for instance in molecu- 
lar quantum dots where stability is a key issue and trans- 
port is the only local probe available. Adiabatic trans- 
port through interacting nano-systems operated in the 
non-linear regime is thus a promising topic where new 
experiments can be done. We acknowledge S. Das, J. 
Konig, M. Plethyukhov, H. Schoeller, C. Stampfer, and 
the financial support from the Helmholtz Foundation, the 
FZ- Jiilich (IFMIT) and the Ministry of Innovation NRW. 
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